c******************************************************************
       function ph(xmu)

       implicit real*8(a-h,o-z)
       ph    =0.5802-0.1571*xmu+.0327*xmu*xmu-.0026*xmu*xmu*xmu
       return
       end

